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Abstract 

Bound state of two massive constituent gluons is studied in the potential approach. Relativistic 
quasi-classical wave equation with the QCD-inspired scalar potential is solved by the quasi-classical 
method in the complex plane. Glueball masses are calculated with the help of the universal mass 
formula. The hadron Regge trajectories are given by the complex non-linear function in the whole 
region of the invariant variable t. The Chew-Frautschi plot of the leading glueball trajectory, ap(t), 
has the properties of the i-channel Pomeron, which is dual to the glueball states in the s channel. 
The imaginary part of the Pomeron is also calculated. 

Pacs: 12.39.Mk; 12.39.Pn; 12.40.Nn; 12.40.Yx 

Keywords: Pomeron, glueball, meson, quark, gluon, potential 

1 Introduction 

Quantum Chromo Dynamics allows the existence of purely gluonic bound states, giueballs. 
These are particles whose valence degrees of freedom are gluons where the gauge field plays 
a more important dynamical role than in the standard hadrons. Study of giueballs is a 
good test of our understanding of the non-perturbative (NP) aspects of QCD, but no firm 
experimental discovery of such gluonic states has been obtained yet. A comprehensive review 
devoted to the giueballs was given in [TH3]. 

An important theoretical achievement in this field has been the computation of the 
glueball masses in lattice QCD [4J. Lattice QCD has been able to compute the low-lying 
glueball spectrum with a good accuracy. This predictions for giueballs are now fairly stable, 
at least when virtual quarks are neglected. Several candidates for the low mass giueballs 
with quantum numbers J PC = ++ , 2 ++ , h and 1 are under discussion. It is difficult to 
single out which states of the hadronic spectrum are giueballs because we lack the necessary 
knowledge to determine their decay properties. 

Many theoretical approaches contain basic model assumptions which are difficult to prove 
starting from the QCD Lagrangian. Complete understanding of giueballs includes such 
theoretical treatments as lattice QCD, constituent models, AdS/QCD methods, and QCD 
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sum rules [J_J[3j|5]. Glueballs have been studied by using effective approach like Coulomb 
gauge QCD and potential model [6HH]. The potential model is very successful to describe 
bound states of quarks. It is also a possible approach to study glueballs (5j[7]. Recent results 
in the physics of glueballs with the aim set on phenomenology and the possibility of finding 
them in conventional hadronic experiments have been reviewed in [5] 

A possible way to handle glueballs is to consider massive quasi-gluons interacting via a 
QCD inspired dynamics. The gluons are massless to all orders in perturbation theory, but 
NP effects like confinement, and their self-interactions, can be described by a constituent 
gluon mass. Another definition of the gluon mass was considered in [9], where a realistic 
QCD motivated gluon propagator was obtained from approximate solution of the Dyson- 
Schwinger equation. The dynamical mass of gluon is defined by the position of the pole of 
the dressed gluon propagator. 

One of open topics in hadron physics is the Pomeron. What is the Pomeron? We know 
the Pomeron as the highest-lying Regge (Pomeranchuk) trajectory (P trajectory) [TD]. In 
the Regge pole theory, the leading Regge trajectories give the main contribution into the 
scattering amplitude. In the many high energy reactions with small four-momentum transfer, 
the soft P exchange gives the dominant contribution in cross sections [IT] . Next question is: 
what is the relation between glueballs and the Pomeron? There has been a long-standing 
speculation that the physical particles on the P trajectory might be glueballs [3j|5j[6]. 

The relation between glueballs and the Pomeron has been investigated by many authors. 
Usually, the P trajectory is considered to be a linear function. However, recent small —t 
ZEUS and HI data for exclusive p and photoproduction [TJHTJ] point out that the P 
trajectory is rather non-linear. The data have been explained by adding in a flavor-blind 
hard Pomeron contribution, whose magnitude is calculated from the data for exclusive J/^f 
photoproduction [T51CEB]. The ZEUS, HI as well as CDF data on pp elastic scattering data 
have also been analyzed by using the non-linear P trajectory [TTH2T)] . 

In this work we consider two-gluon glueballs as the excited states of purely gluonic bound 
states of massive gluons. We accept the potential model, which is so successful to describe 
bound states of quarks; it is also a possible approach to study glueballs [3[2T]. To describe 
the two-body system we use the relativistic quasi-classical (QC) wave equation with the 
QCD-inspired scalar potential, in which the strong coupling is the coordinate dependent, 
i.e., a s = ac s (r). We obtain two exact solutions of the relativistic QC wave equation for 
two components of the potential, the short-distance coulombic term and long-distance linear 
one, separately. Using the interpolation procedure, we join these two solutions and obtain 
an interpolating mass formula, for the bound system. Using this universal mass formula, 
which is good to describe the mass spectra of both light and heavy quarkonia, we calculate 
glueball masses and reconstruct the saturating P trajectory. 

This work is not a comprehensive investigation of the glueball spectroscopy. We concen- 
trate ourself on the leading S z = 2 gluonium states and take a picture where the t-channel 
Pomeron is dual to glueballs in the s channel. We obtain an analytic expression for the P 
trajectory, otp(t), in the whole region of the Mandelstam invariant variable t. The trajectory 
is a complex non-linear function, the real part of which corresponds to the soft Pomeron 
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in agreement with the recent HERA data and is the saturating trajectory. The imaginary 
component of the P trajectory is also calculated. 

2 Glueballs and the Pomeron 

Glueballs were suggested theoretically in [22] and then have been extensively studied in the 
framework of different approaches [21[5j[2Tl[23] • These objects have not been an easy subject 
to study due to the lack of phenomeno logical support. Much debates have been associated 
with their properties. The main achievement of these debates is the understanding of the 
deep relation between the properties of the glueball states and the structure of the QCD 
vacuum. The basic idea is that the vacuum is filled with J PC = ++ transverse electric 
glueballs which form a negative energy condensate [24] . 

Glueballs are bosons made only from the gluonic field; these are quarkless hadrons. They 
can be classed as mesons, because they are hadrons and carry zero baryon number. Glueballs 
must be flavor singlets, i.e., have vanishing isospin (/ = 0) and strangeness. Like all particle 
states, they are specified by the quantum numbers which label representations of the Poincare 
symmetry, i.e., J PC and by the mass. They have the same quantum numbers as isospin 
mesons and their decays in conventional hadrons violate the Okubo-Zweig-Iizuka rule. 

Typically, every quark model meson comes in SU(3) flavor nonets — an octet and a flavor 
singlet. A glueball shows up as an extra (supernumerary) particle outside the nonet. In spite 
of such seemingly simple counting, the assignment of any given state as a glueball remains 
tentative even today. In a strongly coupled theory there is nothing to stop them mixing 
with the quark-based states. These quarkless states are extremely difficult to identify in 
particle accelerators, because they mix with ordinary meson states. Spectrum of pure SU(3) 
Yang-Mills states has been extracted using computerized lattice calculations [I]. 

Pure gauge QCD has been investigated by lattice QCD for many years. This led to a well 
established glueball spectrum below 4 GeV [I]. Lattice QCD has been able to compute the 
low-lying glueball spectrum with a good accuracy. The data shows five isoscalar resonances 
/o(600), /o(980), / (1370), / (1500) and / (1710). Of these the / (600) is usually identified 
with the a of chiral models. The decays and production of /o(1710) give strong evidence 
that it is also a quarkless meson. 

At the present time, the three states 2 ++ , ++ , and h are chosen as possible exper- 
imental glueball candidates [25] because they are computed with relatively small errors in 
lattice calculations [26J. Some experimental glueball candidates are currently known, such 
as the /o(980), /q(1500) and /o(1710), but no definitive conclusions can be drawn concerning 
the nature of these states. We have got no model-independent theoretical knowledge of these 
hadrons. Major new experimental effort forthcoming at Jefferson Lab. 

The modern development in glueball spectroscopy from various perspectives has been 
discussed in [HEKJOj. All these investigations support the use of an effective gluon mass 
to describe the glueball dynamics of QCD. If a valence gluon is a priori massive, then it is 
a spin-1 particle. Such two gluons with total spin S = 2 in glueballs reproduce properly 
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the lattice QCD spectrum for C = + states [I] and may have five 2S + 1 spin states: 
S z = +2, +1, 0, — 1, — 2. Highest glueball trajectory with S z = +2 is expected to be the 
Pomeron. There are works, arguing that a valence gluon is a massless particle, which gains 
a constituent mass fi g , either constant, or state-dependent [7H9]. 

There are other definitions of the gluon mass [IH3], which we do not discuss here. All 
these arguments support the use of an effective gluon mass to describe the dynamics of QCD. 
It is therefore possible to envisage an approach to bound states made of constituent massive 
gluons. We consider here the simplest case of two-gluon glueballs, since they have always a 
positive conjugation charge. 

An open topic in hadron physics is the relation between glueballs and the Pomeron. In 
gauge theories with string-theoretical dual descriptions, the Pomeron emerges unambigu- 
ously. In the QCD framework the Pomeron can be understood as the exchange of at least 
two gluons in a color singlet state [27J. The pQCD approach to the Pomeron, the Balitskii- 
Fadin-Kuraev-Lipatov (BFKL) Pomeron, has been discussed in [2H]. The Pomeron can also 
be associated with a reggeized massive graviton [2"9] . 

The Pomeron is the vacuum exchange contribution to scattering at high energies at 
leading order in 1/N C expansion. It is the highest-lying Regge trajectory. In the many high 
energy reactions with small four-momentum transfer the P exchange gives the dominant 
contribution [IT]. The classic soft Pomeron is constructed from multi-peripheral hadronic 
exchanges. It is usually believed that the soft P trajectory is a linear function, 

a P (t) = a P (0) + ap(0)t, (1) 

where the intercept cnp(O) = 1 and the slope a' P (0) = 0.25 (GeV/c)~ 2 . These fundamental 
parameters are very important in high-energy hadron physics. Usually, they are determined 
from experiment in hadron-hadron collisions. 

To explain the rising hadronic cross sections at high energies, the classic soft Pomeron was 
replaced by a soft supercritical Pomeron with an intercept ap(0) ~ 1.08. The approximate 
linearity ([T]) is true only in a small —t region. The ZEUS, HI as well as CDF data on 
pp elastic scattering data have also been analyzed by using the non-linear P trajectory 
|17j . Important theoretical results have been obtained in [30]. The results imply that the 
effective P trajectory flattens for — t > 1 (GeV/c) 2 that is evidence for the onset of the 
perturbative 2-gluon Pomeron. These results may shed some light on the self-consistency of 
recent measurements of hard-diffractive jet production cross sections in the UA8, CDF and 
HERA experiments. 

The issue of soft and hard Pomerons has been discussed extensively in the literature 
[T0~ | l2"8" | l3~l~] . Both the IR (soft) Pomeron and the UV (BFKL) Pomeron are dealt in a unified 
single step. On the basis of gauge/string duality, the authors describe simultaneously both 
the BFKL regime and the classic Regge regime ■ The problem was reduced to finding 

the spectrum of a single j-plane Schrodinger operator. The results agreed with expectations 
for the BFKL Pomeron at negative t, and with the expected glueball spectrum at positive 
t, but provide a framework in which they are unified. 
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A model for the Pomeron has been put forward by Landshoff and Nachtmann where it is 
evidenced the importance of the QCD NP vacuum [10]. The current data is compatible with 
a smooth transition from a soft to a hard Pomeron contribution which can account for the 
rise of Otot with s. If soft and BFKL Pomeron have a common origin, the discontinuity across 
the cut in the cep(t) plane must have a strong t dependence that points out on non-linearity 
of the P trajectory [32J. 

On the theoretical front, Tang |33j used perturbative QCD (pQCD) to show that Regge 
trajectories are non-linear by studying high energy elastic scattering with mesonic exchange 
in the case of both fixed and running coupling constants. On the experimental side, Brandt 
et al. [31] affirmed the existence of non-linear P trajectories from the data analysis of the 
UA8 and ISR experiments at CERN. They published a parametrization of P trajectories 
containing a quadratic term, 

a P (t) = 1.10 + 0.25* + a' P (0)t 2 , (2) 

where a P (0) is a constant, which are found from the Pomeron data fit. 

Burakovsky et al. [35] presented a phenomenological string model for logarithmic and 
square root Regge trajectories. They applied a phenomenological approach based on non- 
linear Regge trajectories to glueball states. The parameters, i.e., intercept and threshold, or 
trajectory termination point beyond which no bound states should exist, were determined 
from Pomeron (scattering) data. They predicted masses of glueballs on the tensor trajec- 
tory. The approach was applied to available quenched lattice data and found a discrepancy 
between the lattice based thresholds and the Pomeron threshold that was extracted from 
data. 

Linear trajectories are, in fact, disfavored by various experimental data. For more details 
see discussions in [35]. The square root form of the trajectory with the parameters fitted to 
scattering data alone gives the same mass predictions as the fit to both the scattering data 
and the tensor glueball mass, but with larger errors. Using the fit, the authors obtained the 
following predictions for excited glueball masses: M(2 ++ ) = 2.38 ± 0.12 GeV, M(4 ++ ) = 
4.21+0.21 GeV, and M(6 ++ ) = 5.41 + 0.28 GeV with the same central value obtained from 
purely scattering Pomeron data. 

Glueballs in full QCD are very complicated systems. The lightest glueballs with positive 
charge conjugation C = + can be successfully modeled by a two-gluon system (gluonium) in 
the framework of the potential approach [5]. Gluonium leading states and their connection 
with the Pomeron have been studied in our ref. |20j. We modeled glueballs to be bound 
states of two constituent massive gluons interacting by the Cornell potential and shown a 
good agreement both with the lattice calculations in bound state region and the scattering 
data for the Pomeron. 

The funnel-shaped Cornell potential is fixed in an extremely simple manner in terms of 
very small number of parameters. In pQCD, as in QED the essential interaction at small 
distances is instantaneous coulombic one-gluon exchange (OGE); in QCD, it is qq, qg, or gg 
Coulomb scattering [32]. Therefore, one expects from OGE a Coulomb-like contribution to 
the potential, i.e., Ks( r ) ^ —o. s / r at r — > 0. 
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For large distances, in order to be able to describe confinement, the potential has to 
rise to infinity. From lattice-gauge-theory computations [36] follows that this rise is an 
approximately linear, i.e., Vi(r) ~ or+const for large r, where a ~ 0.15 GeV 2 is the string 
tension. These two contributions by simple summation lead to the famous funnel-type 
(Cornell) quark- ant iquark potential [36H38] . 



4 a 

V(r) = V s (r) + V L (r) = --— + or; 

3 r 



(3) 



its parameters are directly related to basic physical quantities noted above. All phenomeno- 
logically acceptable QCD-inspired potentials are only variations around this potential. 

As for gluonium, the situation is very similar. The potential of gg interaction has a 
similar form, but different parameters. A new method called the Vacuum Correlator Model 
(VCM) has been used in [8]. In this model all NP and perturbative dynamics of quarks and 
gluons is universally described by lowest cumulants, i.e., gauge invariant correlators of the 
type (F^xx) . . .F\ a (x u )). In the adjoined and fundamental representations, the final form 
of interaction of two massive gluons is the funnel- type potential of the form fl3]) [8]: 



n/ 

V(r) = -^ + o a r + C , 
r 



(4) 



where a a 
coupling, a 



a 

fund 



adj 



3a{. und , a a = a ad i = (9/4) a^ und ; here a{ und is the light quarks strong 
o — 0.15 GeV 2 is the string tension. The potential @ was used in [20] to 



calculate glueball masses and the P trajectory. 

It is hard to find an analytic solution of the wave equation for the potential (HI). How- 
ever, joining two exact solutions obtained separately for the short-distance coulombic, Vs(r), 
and long-distance linear part, Vi(r), of the potential, with the help of the two-point Pade 
approximant we obtained the interpolating mass formula [20, 39"]|4"0]. 



Mi = 8a 2n 



J + --5 



+ Am 



a 2 



(n r + J + 1 



(5) 



where parameters depend on the system (qq or gg): a = (4/3)af und (qq), 3a{. und (gg) and 
a = <jf und (qq), (9/4)a^ und (gg). The simple mass formula (JSJ describes equally well the mass 
spectra of all qq and QQ mesons ranging from the ud (dd, uu, ss) states up to the heaviest 
known bb systems [3H1SD] ■ This same formula has been used to calculate the glueball masses 
as well [H] . 

Regge trajectories are usually assumed to be linear in t, but there are both phenomenolog- 
ical and theoretical arguments supporting the idea of non-linear trajectories [35]. Inverting 
05]), we obtained the cubic equation for the angular momentum J and, therefore, the analytic 
dependence J(M 2 ) for Regge trajectories including the P trajectory, acp(t), in the whole re- 
gion of the invariant variable t [2011311110] • The "saturating" Regge trajectories in fig. [TJ were 
applied with success to the photoproduction of vector mesons that provide an excellent si- 
multaneous description of the high and low —t behavior of the jp — > pp, u, <p cross sections, 
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Figure 1: The Chew-Frautschi plots of the leading p, eft and P complex Regge trajectories calculated 
from the mass formula ([5]) with the parameters found from the fit of combined HERA p and <j) data 
(triangles) [T2HT3] . and 2 ++ glueball candidate /o(1710) (cross) [IT] . Solid curve is the P trajectory 
given by (|5|), dotted and dashed lines show the leading p and (j) saturating Regge trajectories. 



given an appropriate choice of the relevant coupling constants (JML-model) jlSHU]. As 
was explained in [I5JII6] the hard-scattering mechanism is incorporated in an effective way 
by using the "saturated" Regge trajectories that are independent of t at large momentum 
transfers [20l[39lIlQ] . 

Saturating trajectories have a close phenomenological connection to the quark- ant iquark 
interaction which governs the mesonic structure [20, ESS HO] . They provide an effective way 
to implement gluon exchange between the quarks forming the exchanged meson [J7, 48J and 
lead to the asymptotic quark counting rules [H] that, model independently, determine the 
energy behavior of the cross section at large —t. This approach was successfully adopted to 
explain the large momentum transfer hadron-hadron interactions, as well as several photon- 
induced reactions [IB]. The pion saturating trajectory (a^ at (t) = —1 when t — > — oo) is in a 
form that reproduces the jp — > mr + reaction around #* = 90 |46j . 

A fair agreement with the experiments is achieved when saturating Regge trajectories [201 
[391110] are used for the propagators of the various exchanged mesons. This is an economical 
way to deal with hard scattering mechanisms since the saturation of the Regge trajectories 
(approaching —1 when —t — > oo) is closely related to the OGE interaction between quarks 
[391 SO]. The to meson production channel is particularly instructive in this respect since 
pion exchange dominates the cross section [H]. This same technics and the mass formula 
([5]) have been used in |20j to reconstruct the saturating P trajectory, ap(t), in the whole 
region of t. Our P trajectory and its comparison with known others are shown in fig.[2l The 
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Figure 2: The vacuum effective J = 1 + 2 Regge trajectory. Solid curve is the P trajectory obtained 
from ([5]) with the parameters found from the fit to combined HERA p (triangles) and (p (circles) 
data |12H14j . and 2 ++ glueball candidate fo(17W) (cross) [H]. Other lines show the classic "soft", 
BFKL, and Donnachie-Landshoff "hard" Pomerons. 



P trajectory in figs. [T] and [2] corresponds to the soft Pomeron and is in agreement with the 
HERA data [FMH] : it is the saturating Regge trajectory, i.e., ap(t — > — oo) — > — 1. 

However, the strong coupling a s in the Cornell potential above and in the mass formula 
05]) is a constant value (free parameter). As we know, the strong coupling in QCD is the 
function of q 2 = t, i.e., a s = a s (q 2 ) is the running strong coupling. Below, we introduce the 
dependence a s (r), and obtain similar mass formula and Regge trajectories with the use of 
the QCD-inspired potential. 

3 The QCD-inspired potential 

The strong coupling in QCD is a function of the squared four- momentum transfer, t: a s (q 2 = 
t), or a s (r) in the coordinate space. A more accurate calculation of hadronic masses and their 
trajectories requires the accounting for the dependence a s (r) in the potential of interaction. 

To find the dependence a s (r) let us consider the concept of dynamically generated gluon 
mass, which arises from an analysis of the gluon Dyson-Schwinger (DS) equations [9J. The 
infinite set of couple DS equations cannot be resolved analytically. One must resort to a 
truncation scheme. Cornwall found a gauge-invariant procedure to deal with these equations 

An approximate resolution of the DS equations was obtained in the Feynman gauge. In 
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Euclidean space, this solution is given by D^ u = —ig^ lu D(q 2 ), where, 

q 2 + /! 2 (<? 2 ) ' 



*oD(q 2 ) = (6) 



and 



a M ) = —T- = l T-TO I /I / A 91 ( 7 ) 



4tt 6oln{[9 2 + V(g 2 )]/A 2 } 
with the momentum dependent dynamical mass given by 



fx 2 (q 2 ) = [i 



ln[(g 2 + 4^)/A 2 
ln(V/A 2 ) 



12 
11 



Here in eqs. (J6])-(IE]) bo = (33 — 2n/)/(127r), n/ is number of flavors, /i s = /i(0), A is the QCD 
dimensional parameter; typical values are: fi g = 500 ± 200 MeV and A = 300 MeV [9]. 

This solution contains a dynamically generated gluon mass (IE]) and is another NP ap- 
proach which has led to a very appealing physical picture establishing that the QCD running 
coupling freezes in the NP regime. Expression (J7J is considered to be the QCD running cou- 
pling in momentum representation; it is frozen in the NP regime (q 2 — > 0), 

a = a s (0) = j- — — -, (9) 

because of the presence of the dynamical gluon mass the strong effective charge, g{q 2 ), 
extracted from these solutions freezes at a finite value, giving rise to an infrared fixed point 
for QCD [9] . The gluon mass generation is a purely NP effect associated with the existence 
of infrared finite solutions for the gluon propagator. 

Solution ([6]) is valid only for \i g > A/2. An important feature of the propagator ([6]) is that 
it incorporates the correct ultraviolet behavior, i.e., asymptotically obeys the renormalization 
group equation. This means that the gluon propagator (jSJ) asymptotically at large q 2 takes 
the usual form, i.e., D(q 2 ) oc 1/g 2 , because m 2 (q 2 ) — > at q 2 — > oo and is valid for the entire 
range of momentum. The gluon is massless at the level of the fundamental QCD Lagrangian, 
and remains massless to all order in pQCD. The NP QCD dynamics generates an effective, 
momentum-dependent mass, without affecting the local SU(3)c invariance, which remains 
intact [H]. 

According to Low and Nussinov the Pomeron is modeled as the exchange of two gluons 
[27]. In [50] we modeled the P exchange by two NP gluons as suggested by Landshoff and 
Nachtmann |51j. We dealt with the Cornwall's propagator in the context of the Landshoff- 
Nachtmann model and extracted the "pure" NP propagator. We have shown that the last 
one in combination with the multi-Pomeron asymptotic of the Quark-Gluon String Model 
(QGSM) [52H54"] . results in a good description of soft and hard distributions of secondary 
hadrons in a wide energy range. 

Some consequences of the Cornwall's solution for the gluon propagator associated to 
the static interaction were investigated in [52]. The OGE static potential derived from the 



9 



DS equations (DS potential) was calculated numerically and compared to phenomenological 
potentials whose shape has been inspired by lattice computations. Application of this DS 
potential and comparison with some others to the description of quarkonia was considered. 

The strong coupling a s in the Cornell potential (jHJ), (BJ, is a free parameter. This potential 
can be modified by introducing the a s (r)-dependence, which is unknown. However, using 
the mnemonic rule, q 2 — > 1/r 2 , from (J7|) one can write an ansatz, for the strong running 
coupling in the coordinate space as follows: 

as{r) = 6 ln[l/(Ar)2 + (2/x fl /A)2]- (10) 

The running coupling (fTUj) conserves the basic properties of the one ([7]) in the momentum 
representation: a s {r — > 0) = (q 2 —> oo) and a s (r — > oo) = «o (q 2 —> 0). We see, that 
the running coupling (flOj) is frozen in the NP regime (r — > oo) and is in agreement with the 
asymptotical freedom properties [a s (r — > 0) — > 0]. 

Thus, with the help of (llOp . we come to the following potential of interaction: 

V(r) = -^ + ar, (11) 
r 

where a{r) = ka s (r), k = 4/3 (qq systems), k = 3 (gg system) and a as in ([5]). The spin- 
dependent corrections to the potential (fTTj) can also be derived from lattice QCD, but we do 
not consider them here. 

In hadron physics, the nature of the potential is very important. There are normalizable 
solutions for scalarlike potentials, but not for vectorlike [56] . No any problems arise and no 
any difficulties encountered with the numerical solution if the confining potential is purely 
scalarlike. The effective interaction has to be Lorentz-scalar in order to confine quarks and 
gluons [56J. 



4 Solution of the QC wave equation 



It is hard to find an analytic solution of known relativistic wave equations for the potential 
f fTT]) . that does not allow us to get an analytic dependence E 2 (n r ,l). This aim can be 
achieved with the use of the QC wave equation [57l[58], which for two bound particles of 
equal masses in the rest frame is 



+ 



if 



+ 



1 d 2 



d6 2 sin 2 6 dp 2 



+ 



£2 
4 



(m + Vf 



i)(f) = o. 



;i2) 



This is the second-order differential equation of the Schrodinger type in canonical form. 
Important feature of this equation is that, for two and more turning-point problems, it can 
be solved exactly by the conventional WKB method [571459] . 
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Appropriate solution method of the QC wave equation, which is the same for relativistic 
and non- relativistic systems, was developed in [58j[60]. In our method, each of the one- 
dimensional equations obtained after separation of the QC wave equation is solved by the 
same QC method. The QC wave equation (fl2i) is separated that gives, for the potential (TTTh . 



dr 2 



El 
4 

' 

de 2 



air) 



m 



+ ar 



M 2 



R{r) = 0, 



+ 



d 2 



sin 9 dip' 



+ M Z 



Y(6,<p) = 0. 



(13) 



(14) 



The angular QC equation (|14p determines the squared angular momentum eigenvalues, M 2 , 
which enter into the radial equation f JT3|) . Solution of ffT4"j) has been obtained in [581159] by the 
QC method in the complex plane that gives M 2 = {1 + \) 2 - This result means that the radial 
QC wave equation has always the centrifugal term {I + \) 2 /r 2 for all spherically symmetrical 
potentials V{r). The squared angular momentum eigenvalues, {I + |) 2 , are universal for all 



central potentials and not any Langer-type corrections are required 

This QC method reproduces the exact energy eigenvalues for all known solvable problems 
in quantum mechanics [58] • In our QC method not only the total energy, but also momentum 
of a particle- wave in bound state is the constant of motion. Solution of the QC wave equation 
in the whole region is written in elementary functions as [60J, 



R{r) 



J_ p \Pn\r- 



r < ri, 



C n { cos(|p n |r - 0i - f), ri < r < r 2 , 



V2 



(15) 



r > r 2 , 



where C n = y2|p n |/[7r(n + ^) + 1] is the normalization coefficient, p n is the corresponding 
eigenmomentum, (pi = —n(n + |)/2 and 02 = 7r(n + |)/2 are the values of the phase-space 
integral at the turning points x\ and X2, respectively. In the classically allowed region [x±, X2], 
the solution is 

Rniij) = C n cos (\p n \r + ^nj , (16) 

i.e., has the form of a standing wave. This solution is appropriate for two-turning-point 
problems both in non-relativistic and relativistic cases with the corresponding eigenmomenta 

Pn- 

For example, in case of the non-relativistic Coulomb problem, the total energy eigenvalues 
have the form of kinetic energy of a free particle [58,60j, 



E„ 



Pn_ 

2m' 



Pn 



n r + l + V 



(17) 



where p n = mv n is the non-relativistic momentum eigenvalue with the imaginary discrete 
velocity, v n = ia/(n r + / + 1). This means, for example, that the motion of the electron in 
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a hydrogen atom is free, but restricted by the "walls" of the potential... This is free finite 
motion of a particle-wave in bound state. One should note, that the QC eigenf unctions 
( JIBjl correspond to the asymptote of the exact solution of the Schrodinger equation, i.e., the 
principal term of the asymptotic series of the corresponding exact solution [60J. 

The radial QC equation ( fl~3l) has four turning points and can not be solved analytically by 
standard methods. Let us use the QC method to solve the equation. The QC quantization 
condition appropriate to (TT51) in the complex plane is 




(/ + ^dr = 4n(n r + -) . ( l> 



To calculate the phase-space integral ([IB]) in the complex plane we chose a contour C en- 
closing the cuts (therefore, turning points and zeros of the w.f.) at r < and r > between 
the turning points 7*1, r 2 and r 3 , r 4 , respectively. Outside the contour C, the problem has 
two singularities, i.e. at r = and oo. Using the standard method of stereographic projec- 
tion, we should exclude the singularities outside the contour C . Excluding these infinities 
we have, for the integral ffT8l) , I = I + J^, where Jo = —2n(l + |) is contribution of the 
centrifugal term. The integral 1^ is calculated with the help of the replacement of variable, 
i.e., z — 1/r, that gives = 2n(E 2 /8a + a ), where a = ka and a is the strong coupling 
in the NP regime OH]). Here we took into account the asymptotic properties of the running 
coupling (flU!) and its derivative: a s (r) = 0, ct' s (r) = at r — > 0. The calculations result in 
the squared total energy eigenvalues, 

E n (J) 2 = 8a(2n r + J + 3/2-a ). (19) 

Putting in (Tl9|) « = we come to well known relativistic result for the linear potential, 
E n (J) 2 = 8a(2n r + J + 3/2). 

It is an experimental fact that the dependence E 2 (J) is linear for light mesons. However, 
at present, the best way to reproduce the experimental masses of light hadrons is to rescale 
the entire spectrum assuming that the masses M n of the mesons are expressed by the relation 

m 

Ml = E 2 n - C 2 , (20) 

where C is a constant energy (free parameter). Relation (12"0"|) is used to shift the spectra 
and appears as a means to simulate the effects of unknown structure approximately. This 
constant can be interpreted as a renormalization of the vacuum energy [62]. It has been 
suggested that the confinement potential has a complex Lorentz structure, and the relation 
(|20l used to shift the spectra appears as a means to simulate approximately the effects of 
this structure. This constant can be connected with the structure of the QCD vacuum, filled 
with J PC = ++ transverse electric glueballs which form a negative energy condensate [24]. 

The oscillator-type expression ffl9|) does not require any additional free parameter. It 
contains the needed shift in the form of the interference term — Sao^ of the coulombic and 
linear terms of the potential (TTT]) - Formula ffl9|) is good to describe spectra of light hadrons, 
but not heavy quarkonia. 



12 



The asymptotic expression (|T9|) is defined by the singularities at r = and r — > oo. 
The leading singularity in the phase-space integral ( ITS]) at r — > oo is given by the quadratic 
term (err) 2 originating from the linear part of the potential ( ITT1) . which gives the dominant 
contribution in ( fl8l) at infinity. But the term oc r 2 suppresses a very important contribution 
of the coulombic interaction at small and moderate distances r. 

It is known that heavy QQ systems can be treated non-relativistically and, for low states 
of heavy quarkonia, the main contribution to the boundary energy comes from OGE term 
of the potential (TIT]) , i.e., in the first approximation one can neglect the confining linear 
term. A reliable consideration of the excited QQ states requires a completely relativistic 
treatment. 

Let us consider the radial equation ( TT3|) just for the coulombic part Vs(r) of the potential 
( ITT]) . The QC quantization condition in the complex plane is 



c 



E 2 
T 



a(r) 



m 



dr = 2tt [n r + 



(21) 



and the integral is calculated analogously to the above case. A contour C encloses the 
classical turning points ri and r 2 and cut between them. Using the same solution method 
of stereographic projection, we obtain, for the integral fl2Tj) outside the contour C at r = 
and oo: I = I + 1^, where I = —2n(l + |) and 1^ = 2ir(a m/ J—E 2 /A + m 2 . This gives: 



El = Am 2 



a 



n r + I + 1 



+ 1 



(22) 



Again, here we took into account the asymptotic properties of the strong running coupling 
(|T0|) at r and r — s- oo. 

Thus, we have two exact analytic expressions (I19p and (122 p for two asymptotic compo- 
nents of the potential (TTTj) (coulombic and linear). Now, we can use the same approach as 
in [20l|39lll0] , i.e., derive the interpolating mass formula for E 2 , which satisfies both of the 
above constraints: the exact energy eigenvalues (TT9|) and fl22|) . To derive such a formula we 
use the two-point Pade approximant [63] . 



[K/N] 



/I 



diZ" 



(23) 



with K = 3 and N = 2. We take K — 3 and N = 2 because this is a simplest choice to 
satisfy the two asymptotic limits ( IT9l) and ( )22l) . Simple calculations give the interpolating 
mass formula, 



Mi 



Am 2 



( 2;/, 



T 3 

J+ 2 



a 



n r + J + 1 



(24) 



Note, the two exact asymptotic expressions f|T9|) and f )22|) for E 2 have the form of the squared 
total energy for two free relativistic particles, i.e., E 2 = Aij)^ + m 2 ). 
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The mass formula ( 124)) is good to describe the mass spectra both light and heavy quarko- 
nia. To demonstrate its efficiency we calculate the leading state masses of p and families 
(see tables 1, 2, where masses are in MeV). 



Table 1. The p-family J = I + 1 leading states 



Meson 


jPC 


rpex 


Ef 


Parameters in (1241) 


p{\S) 


1— 


775 


775 


A = 487 MeV 


a 2 {lP) 


2++ 


1318 


1319 


a = 0.137 GeV 2 


Ps(lD) 


3~ 


1689 


1689 


m n = 144 MeV 


a A (lF) 


4 ++ 


2001 


1989 




P (1G) 


5— " 




2249 




P (1H) 


Q++ 




2481 




p(2S) 


\ — 


1720 


1683 




p(2P) 


2++ 




1986 




p(2D) 


3~ 




2247 




p(2F) 


4 ++ 




2480 




p(3S) 


1~ 




2245 




Table 2. The 


0-family J = I 


+ 1 leading states 


Meson 


jPO 


Tpex 


E* 


Parameters in (1241) 


<I>{1S) 


1™ 


1020 


1019 


A = 445 MeV 


Mip) 


2++ 


1525 


1525 


a = 0.125 GeV 2 


Mm 


3" 


1854 


1854 


m s = 414 MeV 




4 ++ 


2018 


2119 




4>(1G) 


5~ 




2349 




<I>{1H) 


6++ 




2556 




4>(2S) 


1~ 


1820 


1820 




4>{2P) 


2++ 


2011 


2103 




<P(2D) 


3~ 




2340 




<j>{2F) 


4 ++ 




2551 




<f>(3S) 


1~ 




2327 





In this calculations the frozen strong coupling (|9]) depends on the ratio p g /A. From the 
fit results we found the optimal value for the constituent gluon mass, p g = 416 MeV. The 
QCD dimensional parameter A as others is found from the best fit to the available particle 
data [64]. For light quarks, we take the average effective mass, m n = (m u + md)/2 MeV. 
The best fit to the data is achieved for m n = 144 MeV. 

The mass formula (1241) is appropriate to calculate the glueball masses as well. In this 
case we use the potential (TTTT) with the parameters of gg interaction. Calculation results for 
the gluonium leading state masses are shown in table 3. 
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Table 3. Glueball J = I + 2 leading states 



Glueball 




rpex 


rpth 


Parameters in (124 jl 


/o(1710) 


2++ 


1710 


1710 


A = 321 MeV 


fo (IP) 


3" 




2405 


cr a = 0.329 GeV 2 


fo (ID) 


4++ 




2921 


ti g = 416 MeV 


fo (IF) 


5~ 




3350 




fo (2S) 


2++ 




2898 




fo (2P) 


3~ 




3338 




fo (2D) 


4++ 




3720 




fo (3S) 


2++ 




3711 




fo (SP) 


3— " 




4057 




fo (3D) 


4 ++ 




4374 





Equation (I24p is an ansatz [as the potential (Hip ], which is based on two exact asymptotic 
expressions ( 1T9|) and (|22|) . It allows us to get an analytic expression for Regge trajectories 
in the whole region. Transform (j24p into the cubic equation for the angular momentum 
J(t = E 2 ), 

J 3 + ai (t)J 2 + a 2 (t)J + a 3 (t) =0, (25) 

where a x (t) = X(t) + 2n, a 2 (t) = 2nX(t) + n 2 , a 3 (t) = n 2 X(t) - a 2 m 2 /(2a), n = n r + 1, 
X(t) = (—t + 4m 2 )/(8a) + 2n r + 3/2 — a. Equation (|25|) has three (complex in general case) 
roots: J\(t), J 2 (t), and J%(t). The real part of the first root, Re J\(t) (Chew-Frautschi plot), 
gives the analytic expression for Regge trajectories, 

f 2 v ^cos(0/3) - ai/3, Q < 0; 
Rea(t) = l -ai/3(g = 0), Q = 0; (26) 

l/i + / 2 -ai/3, g>0, 

where 0(t) = arccos[— qf \J— p' 6 ], p(t) = — a 2 /9 + a 2 /3, q(t) = af/27 — axa 2 /6 + a 3 /2, Q(t) = 
P 3 + Q 2 , fi(t) = [-Q + v / Q] 1//3 ) = [ — Q — VQ] 1 ^ 3 - Expression f[26]) supports the existing 
experimental data and gives the saturating Regge trajectories including the Pomeron in the 
whole region of t. 

The imaginary part of the complex Regge trajectories is given by the expression 

lma(t) = V3(f 1 -f 2 )/2, Q>0. (27) 

The threshold (trajectory termination point), beyond which no bound states should exist, is 
defined from the equation, Q(t) = p 3 + q 2 = (see notations above). 

The effective P trajectory in figs. [U13] has similar properties as all quark- ant iquark tra- 
jectories given by ( )26l) . It is asymptotically linear at t — > 00 with the slope a! P = l/(8er a ) ~ 
0.380 (GeV/c) -2 , and flattens off at —1 for t — > —00. As all Regge trajectories the P trajec- 
tory is the monotonically rising function at the interval (—00, 00). 

Parameters of the function Reap(t) are found from the best fit to the combined ZEUS 
p (triangles) and (circles) scattering data [12], and glueball candidate /o(1710) [H] (cross) 
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with quantum numbers J PC = 2 ++ (/ = 0, S z = +2) in bound state region. The intercept 
and slope of the P trajectory at t = are: 



a'p(O) = 1.083, a'p(0) = 0.280(GeV/c) 2 . 



(28) 



The hard BFKL Pomeron has intercept «p(0) ~ 1.43 [28]. The Pomeron with such properties 
results in too fast growth of the total cross sections. Recent studies of exclusive electro- 
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Figure 3: The real and imaginary components of the gluonium J = 1 + 2 complex Regge trajectory 
(the Pomeron). Triangles show combined HERA p and (j) data |12H14j . 



production of vector mesons at JLab [121113] made it possible for the first time to play with 
two independent hard scales: the virtuality Q 2 of the photon, which sets the observation 
scale, and the momentum transfer t to the hadronic system, which sets the interaction scale. 
They reinforce the description of hard scattering processes in terms of few effective degrees 
of freedom relevant to the Jlab- Hermes energy range [47] . 

The study of exclusive electro-production of u> mesons, completed at JLab [12JH3], pro- 
vides us with an original insight on the space time structure of hard scattering processes 
between the constituents of hadrons. It was shown that the higher order mechanisms are 
more economically described in terms of a few effective degrees of freedom: dressed par- 
ton propagators, saturating Regge trajectories and electromagnetic form factors of off-shell 
meson. The success of this description in several channels is a strong hint that they are 
the relevant degrees of freedom in the JLab-Hermes energy range. In addition, they pro- 
vide us with a link with more fundamental approaches of NP QCD: ab initio Lattice Gauge 
calculations or potential models. 

The two-Pomeron picture (soft plus BFKL Pomeron) gives a very good fit to the total 
cross section for elastic J/\I/ photoproduction and the charm structure function F% over the 
whole range of Q 2 = —t [65]. The hard BFKL Pomeron has the intercept «bfkl(0) — 1.43. 
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Next-to-leading order estimates give, for the BFKL intercept values 1.26 to 1.30, which is 
closer to the soft supercritical Pomeron. However, the results of the experiments and the 
found higher order corrections make it quite unclear what the BFKL Pomeron is. Another 
question is: what is the intercept of the BFKL Pomeron, if the pQCD is non-applicable 
in this scale, i.e., at t — 0? On the other hand, the saturated P trajectory ( 1261) has the 
properties of the soft supercritical Pomeron at small —t and saturates at large —t according 
to the pQCD prediction. 

5 Conclusion 

We have considered glueballs as bound states of constituent massive gluons and investigated 
their properties in the framework of the potential approach. The constituent gluon picture 
could be questioned since potential models have serious difficulties in reproducing all the 
currently known lattice QCD data. In spite of non-relativistic phenomenological nature, the 
potential approach can be used to describe glueballs. Fair description of quarkonium states 
give us a confidence that we are on the right way to describe glueballs. 

The physical properties of constituent gluons are still a matter of controversy. Within 
the framework of potential models, gluons are supposed to be massless or massive, i.e., with 
either a helicity-1 or a spin-1 particles. We have dealt with simplest two-gluon glueballs, 
but mass and spin in different works are very different. If valence gluons are assumed to be 
helicity-1 particles, then their spin has only two projections: S z = ±1. In this work gluons 
are considered as massive spin-1 particles with projections S z = —1,0,-1-1. The properties 
of these pure glue states are not completely understood. 

We have analyzed two exact asymptotic solutions of relativistic QC wave equation for the 
QCD-inspired scalar potential with the coordinate dependent strong coupling, a s (r). One 
needs to stress that the behavior of the potential and exact form of the coupling a s (r) are 
not so important in the intermediate region: in fact, the results depend on the asymptotic 
behavior of the potential at r — > and r — > oo. Our QC method in the complex plane 
to resolve the eigenvalue problem is rather simple and allows to get the analytic result. 
Using two asymptotes, corresponding to the short-distance coulombic and long-distance 
linear components of the potential, we have derived the universal mass formula and 
calculated the glueball masses, which are in agreement with the lattice data. 

We have considered glueballs as the physical particles on the P trajectory. To reproduce 
the trajectory, we have inverted the mass formula ( |24l) and derived the analytic expression 
( 12"B"|) for the P trajectory, Reap(t), and its imaginary part Im«p(i) in the whole region of 
t. In the scattering region, at —t ^> Aq C d, the trajectory flattens off at —1, i.e., it has 
asymptote ap(t — > — oo) = —1 (saturates). In bound state region, at large timelike t, the P 
trajectory is linear in accordance with the string model. 

It is known, that the fixed-number of particles with a potential description can not be 
used for strict relativistic description. Strict description of the Pomeron presupposes multi- 
particle description of the system. For perturbative regime with the Pomeron scattering, the 



17 



dominant contribution comes from the BFKL Pomeron. However, experimental data and 
our simple calculations in the framework of the potential approach support the conception of 
the soft supercritical Pomeron as observed at the presently available energies. The saturat- 
ing Regge trajectories (l26i) obtained in this work effectively include short- and long-distance 
dynamics of constituents. 

In this paper we have not considered helicities and spin properties of gluons. This topic 
has been discussed in details somewhere else [5]. The existing data and simple analysis 
performed in this work confirm the existence of the Pomeron which is the complex non- 
linear function with the properties of soft supercritical Pomeron at small —t and saturates 
Pomeron at large squared momentum transfers. 
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